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Radial Distribution of the Center of Gravity of Random 

Points on a Unit Circle 

F. Scheid 2 

This paper describes some Monte Carlo computations carried out on the Standards 
Electronic Automatic Computer (SEAC) which are related to the problem of random walks. 



In 1905 Pearson [1] 3 suggested the problem of a 
random walk involving n steps of equal length and 
arbitrary direction in the plane. In 1906 Kluyver [2] 
obtained the solution 



i 



P(r,ri)=r J [J (*)] w e/i(rz)da; 

giving the probability of being no farther than r from 
the starting point after n steps of unit length. This 
integral has been calculated by Greenwood and 
Durand [3] for ^=6(1)24, and r=0.5(.5)/i, using 
punched-card tables of the Bessel functions. These 
tables were not extensive enough to handle the values 
71=3,4, and 5. For n — 2 one finds by a simple argu- 
ment, P(r,2) = (2/ir) arcsin (r/2). 

The problem described in the title is essentially 
the same as Pearson's. Its solution for r&=3.4, and 5 
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has been obtained to roughly three decimal places 
by a direct sampling procedure, using SEAC. A set 
of n pseudorandom numbers is generated by using the 
multiplication method of Taussky and Todd [4]. 
A linear transformation converts these into a sample 
from the uniform distribution over ( — ir,ir). The 
center of gravity is found and classified in one of the 
intervals j/32<fl<(j+l)/32, j=0, . . .,31. A 
new sample is then taken. Both frequency table 
and cumulative distribution are printed. Computa- 
tion of Kluyver's integral for these values of n is 
possible but lengthy. If the integrand is machine 
computed as needed, many hours of computer time 
are consumed because of the two independent 
parameters r and n. For limited accuracy the 
sampling approach seems better. The distributions 
for 7i=2 and 6 have also been obtained for accuracy 
checks. Using a theorem of Kolmogoroff [5], a brief 
calculation indicates that, for n=2, the probability 
of a maximum error less than or equal to 0.01 after 
6,000 samples have been taken is approximately 
one-half. A glance at table 1 shows that in fact the 



Table 1. 



32/2 




71 = 2 




n 


=3 


n- 


= 4 


n = 


= 5 


71 = 6 




Freq 


Cum 


Exact 


Freq 


Cum 


Freq 


Cum 


Freq 


Cum 


Cum G-D 


1 


121 


.0107 


.0199 


7 


. 001 


36 


.005 


28 


.004 


.0000 .0000 


2 


133 


.0413 


. 0398 


37 


.007 


87 


.018 


80 


.016 




3 


120 


.0018 


. 0598 


58 


.017 


128 


.038 


158 


.040 




4 


124 


. 0820 


. 0798 


07 


.028 


169 


.063 


185 


.068 




5 


129 


L030 


. 0999 


95 


.043 


209 


.094 


210 


.099 




6 


111 


. 1211 


. 1201 


113 


.061 


192 


.123 


308 


. 140 




7 


123 


. 1411 


. 1404 


141 


.084 


266 


.163 


381 


.203 




8 


115 


.1598 


.1609 


172 


.112 


289 


.207 


361 


.257 


. 2910 . 2932 


9 


129 


.1808 


.1816 


224 


.149 


238 


.242 


382 


.314 




10 


142 


. 2039 


.2023 


336 


.203 


316 


.290 


350 


.367 




11 


123 


. 2240 


. 2234 


466 


.279 


335 


.340 


361 


.421 




12 


138 


.2464 


.2447 


344 


.335 


360 


.394 


384 


.479 




13 


126 


. 2669 


. 2663 


291 


.383 


357 


.448 


358 


.533 




14 


157 


.2925 


.2883 


285 


.429 


365 


.503 


384 


.590 




15 


126 


.3130 


. 3106 


269 


.473 


365 


.558 


374 


.647 




16 


125 


.3333 


.3333 


255 


.514 


405 


.618 


330 


.696 


. 7665 . 7672 


17 


150 


.3577 


.3565 


223 


.551 


353 


.672 


317 


.744 




18 


158 


.3835 


.3803 


189 


.581 


255 


.710 


324 


.802 




19 


135 


.4054 


.4047 


208 


.615 


275 


.751 


273 


.834 




20 


148 


. 4295 


.4298 


185 


.645 


262 


.790 


224 


.867 




21 


157 


. 4551 


.4558 


215 


.680 


182 


.818 


179 


.894 




22 


158 


.4808 


.4826 


197 


.712 


159 


.842 


143 


.916 




23 


173 


. 5(100 


.5106 


183 


.742 


163 


.866 


131 


.935 




24 


190 


.5399 


.5399 


201 


.775 


168 


.892 


109 


.952 


. 9749 . 9732 


25 


191 


.5710 


.5708 


188 


.805 


167 


.917 


93 


.966 




26 


211 


.6053 


. 6038 


183 


.835 


131 


.936 


70 


.967 




27 


197 


. 6374 


. 6393 


163 


.862 


102 


.952 


45 


.983 




28 


247 


. 0770 


.6783 


176 


.890 


87 


.965 


39 


.989 




29 


202 


. 7202 


.7221 


170 


.918 


87 


.978 


36 


.994 




30 


308 


.7703 


.7737 


162 


.944 


76 


.989 


24 


.998 




31 


424 


siv.n 


.8407 


163 


.971 


45 


.996 


12 


.999 




32 


987 


1.0000 


1.0000 


178 


1.000 


27 


1.000 


3 


1.000 


1.0000 1.0000 


Number of samples . . 


6144 


6144 


6656 


6656 




6656 



307 



Table 2. 



32 R 


Freq Cum 


71 = 4 

Freq Cum 


w = 5 
Freq Cum 


1 


4 . 001 
10 . 003 
17 . 006 
19 . 010 
48 . 019 
45 . 028 
76 . 043 
92 . 061 

107 . 082 
160 . 113 
289 . 169 
203 . 209 
203 . 249 

176 . 283 

177 . 318 
185 . 354 

170 . 387 
172 . 421 
176 . 455 
196 . 493 
185 . 529 
196 . 568 
196 . 606 
172 . 639 

184 . 675 

224 . 719 
220 . 762 
236 . 808 

225 . 852 
239 . 899 
239 . 946 
279 1. 000 


9 .002 
29 . 007 
50 . 017 
59 . 029 
71 . 043 
88 . 060 
69 . 073 
109 . 095 

133 . 121 
126 . 145 

151 . 175 
191 . 212 
221 . 255 

231 . 300 

232 . 346 
287 . 402 

248 . 450 
237 . 496 
235 . 542 
235 . 588 
215 . .630 
191 . 667 
229 . 712 
210 . 753 

205 . 793 
183 . 829 
188 . 866 
183 . 901 
167 . 934 

152 . 964 
120 . 987 

66 1. 000 


7 .002 
17 . 006 
21 .011 
40 . 021 
61 . 036 
72 . 053 
99 . 077 
113 . 105 

106 , .131 
120 . 160 
142 . 195 
149 . 231 

177 . 274 

167 .215 
227 . 371 
217 . 424 

242 . 483 
234 . 540 
237 . 598 

207 . 648 

208 . 699 

178 . 742 

179 . 786 

168 . 827 

148 . 863 
156 . 901 
117 .930 
97 . 954 
77 . 972 
65 . 988 
34 . 997 
14 1. 000 


2 


3 


4 


5 


6 


7 _ - 


8 


9. 


10 


11 


12 


13 


14 


15_ 


16 


17_ 


18__. 


19 


20_ 


21 


22 


23_ 


24 _ 


25 . 


26 . 


27 _. 


28 


29 _. 


30_ -_ 


31_ 


32_ 




Number of samples- _ 


5120 


5120 


4096 



Washington, September 16, 1957. 



maximum error is 0.004 for this distribution. Since 
the inverse square root law is operative, to obtain 
four decimal-place accuracy would require perhaps 
two-thousand times as manv samples, or a year of 
SEAC time, 

Parallel results are given in table 2 for samples 
drawn from a triangular distribution over (— ir,ir). 
The pseudorandom numbers are modified by a direct 
method [6] to produce the required sample. 

The question how best to compute a distribution 
function is not, of course, answered in this paper. 
It is suggested, however, that even in cases where 
the solution can be written explicitly, circumstances 
may make computation by direct sampling 
preferable. 
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